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ABSTRACT 

We present an analysis of XMM-Newton and RXTE data from three observations of the neutron star 
LMXB 4U 1636-536. The X-ray spectra show clear evidence of a broad, asymmetric iron emission line 
extending over the energy range 4-9 keV. The line profile is consistent with relativistically broadened 
Fe K-a emission from the inner accretion disk. The Fe K-a line in 4U 1636-536 is considerably broader 
than the asymmetric iron lines recently found in other neutron star LMXBs, which indicates a high 
disk inclination. We find evidence that the broad iron line feature is a combination of several K-a 
lines from iron in different ionization states. 

Subject headings: accretion, accretion disks — binaries: close — stars: individual (4U 1636-536) — 
stars: neutron — X-rays: binaries — X-rays: stars 



1. INTRODUCTION 

Relativistically broadened, asymmetric Fe K-a lines 
from the inner accretion disk have been observed in 
many superm assive and stellar-mass black holes (e.g. 
iFabianI l2006f ). In neutron star binaries, however, 
the iron lines are weaker, and until recently observa- 
tions did not clearly reveal a relativistic line profile. 
[Bhattacliarvva & Strohmayer (2007) found an asymmet- 
ric Fe K-a line in XMM-Newton spectra of the low-mass 
X-ray binary (LMXB) Serpens X-1 and showed that the 
line profile is consistent with fluorescent iron line emis- 
sion from the inner accretion disk. Similar asymmet- 
ric Fe K-a lines were found bv ICackett et al . (2008) in 
Suzaku spectra of the neutron star LMXBs Serpens X-1, 
4U 1820-30, and GX 3A9+2. In this paper we present 
XMM-Newton and RXTE observations of the neutron 
star LMXB 4U 1636-536. We analyze the X-ray spectra 
to determine the profile of the relativistic Fe K-a line 
and constrain the properties of the inner accretion disk. 

4U 1636-536 (V801 Ara) is a well-studied, burst- 
ing LMXB consisting of a neutron star in a 3.8- 
hr orbit with a 0.4 sola r mass, 18th magnitude star 
(|van Paradiis et aLlll990D and is l ocated at a distance 
of ~6 kpc (|Gallowav et all |2006[ ). The X-ray tim- 



ing pro perties of th e bina ry ha ve been studied exten- 
sively. IZhang et al] (|1996l ) and IWijnands efal] (|1997l ) 
discovered quasi-periodic oscillations (QPOs) at kHz 
frequencies. The source also exhibits highly coherent 
burst oscillations at 581 Hz which are likely related 
to the rotation of the neutron sta r (IZhang et a,l1 119971 : 
IStrohmaver fc Markward^ I2002D . iKaaret et aLr(|1999f ) 
found that the soft X-ray emission, modulated at the kHz 
QPO frequency, lags behind the hard X-ray emission. A 
possible explanation for this phase lag is the reprocess- 
ing of hard X-rays in a cooler Comp tonizing corona with 
a size of at most a few kilometers. iBarret et al.l ()2007l ) 
interpreted observations showing a decline of the QPO 
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coherence and rms amplitude at high QPO frequencies 
as evidence that the inner radius of the accretion disk 
in 4U 1636-536 is usually larger than but sometimes 
approach es the iii nermost stable circular orbit (ISCO). 
However, iMendej ([2006) argued that this decline is not 
caused by effects related to the ISCO. 

In this paper we report on three X-ray observations of 
4U 1636-536 that were carried out simultaneously with 
XMM-Newton and RXTE. We present an analysis of the 
X-ray spectrum over the 0.5-100 keV energy range and 
our results from modeling the continuum and relativis- 
tic Fe K-a line emission. We use the Fe K-a line profile 
to derive constraints on the disk inclination and inner 
disk radius. We show that the line profile is not consis- 
tent with a single relativistic Fe K-a line from a neutron 
star accretion disk and that at least two lines from dif- 
ferent ionization states of iron are needed to adequately 
describe the line profile. Finally, we discuss the impli- 
cations of our findings for measurements of neutron star 
radii based on Fe K-a line profiles. 

2. OBSERVATIONS 

4U 1636-5 36 was observed with XMM-Newton 

(iJansen et al.l l2001l ) on 2005 August 29 (Obs. ID 
0303250201), on 2007 September 28 (Obs. ID 
0500350301), and on 2008 February 27 (Obs. ID 
05003 50401). The EPIC PN camera (IStriider et all 
|2001| ) collected 30.0 ks of data starting at 18:24 UT 
during the 2005 observation (hereafter observation 1), 
30.6 ks of data starting at 15:45 UT during the 2007 
observation (observation 2), and 38.6 ks of data starting 
at 04:16 UT during the 2008 observation (observation 
3). The EPIC PN was operated in timing mode and 
with the medium blocking filter. We processed the PN 
data from observations 1 and 2 with the XMM-Newton 
SAS version 7.1.0 using the latest calibration files. The 
data from observation 3 required special processing by 
the XMM-Newton Science Operations Center because 
of the very large number of events. We extracted 
source photons from the timing mode data using CCD 
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rows RAWX = 30-46 (29-45 for observation 3) and 
background photons using RAWX = 2-18. The full 
range of CCD columns RAWY was used. We selected 
only events with PATTERN < 4 (singles and doubles) 
and FLAG = and restricted our analysis to the energy 
range 0.5-12 keV. The average PN source count rate in 
this energy range, excluding X-ray bursts, was 230 s~^ 
for observation 1, 470 for observation 2, and 620 
for observation 3. These rates are below the 800 s"'^ tim- 
ing mode pile-up limit above which the spectral response 
would be deteriorated by photon pile-up {XMM-Newton 
User's Handbook). The background count rate was less 
than 2% of the source count rate and did not exhibit 
significant flaring. We find a strong line-like feature in 
the background spectrum near 0.45 keV. We interpret 
this feature as the electronic noise peak for double 
events that is also present for PN imaging modes but is 
shifted to higher energies in timing mode. 

4U 1636-536 was observed with RXTE (jBradt et all 
119931 ) simultaneously with each XMM-Newton observa- 
tion. The RXTE observations where somewhat longer 
than the XMM-Newton observations and completely 
overlapped with the EPIC PN exposures. No earth oc- 
cultation occurred during these observations and the tar- 
get was observed continuously with RXTE. In this paper 
we consi der only the data obtained with the HEXTE 
detector (jRothschild et all [199811 . We did not use the 
PCA data because of their lower energy resolution com- 
pared to the EPIC PN data and because the PCA spec- 
trum in the 3-12 keV energy range showed an up to 30% 
higher flux than the EPIC PN spectrum. Similar ex- 
cesses of the PCA flux compared to other instruments 
have previously been re ported fe.g. iTomsick et aI1ll999l : 
iCourvoisier et al.ll2003D . It was not possible to correct 
for the higher PCA flux by including a multiplicative 
cross-calibration factor because the flux excess is strongly 
energy dependent (^30% at 3 keV and ^10-15% at 
10 keV). We included in our analysis only those HEXTE 
data that were taken simultaneously to the EPIC PN 
data. The HEXTE data cover the entire duration of the 
EPIC PN exposures with the exception of a single <1 ks 
data gap in each observation. The HEXTE background 
rate and the telescope pointing were stable during both 
observations. We extracted HEXTE data for the energy 
range 20-100 keV using the REX data analysis script. 
For observation 1, we included in our analysis data from 
both HEXTE clusters. For observations 2 and 3, we only 
used the data from cluster B because cluster A was no 
longer able to obtain background measurements. 

3. SPECTRAL ANALYSIS 

For the spectral analysis, we removed all X-ray bursts 
from the data (three for observation 1 and one each for 
observations 2 and 3) . The EPIC PN spectra were binned 
at about 1 /3 of the FWHM detector resolution and were 
fitted simultaneously with the HEXTE data. We did 
not include a multiplicative cross-calibration factor be- 
tween the two instruments. The spe ctral fitti ng was per- 
formed using XSPEC version 12.4 (jArnaudI [19961. For 
the EPIC PN spectrum from observation 3, we find sig- 
nificant residuals between data and model at ~2.3 keV 
near the instrumental Au edge. We attribute these resid- 
uals to a small offset in the photon energy determination 
caused be a slightly inaccurate charge transfer efficiency 



(CTE ) calibration in timing mode (see e.g. ISala et all 
I2007D . We corrected for this energy offset by introducing 
a multiplicative gain correction factor of 1.0020 which we 
determined from a fit. This gain factor, which was fixed 
in our further spectral analysis, completely removed the 
residuals near the instrumental Au edge. No gain cor- 
rection was necessary for the spectra from observations 
1 and 2. 

Initial fitting of the spectra from all three observations 
showed significant residuals between data and model be- 
low 2 keV near the K-shell absorption edges of O, Ne, 
Mg, and Si. We attribute these residuals to an overly 
simplified description of the interstellar absorption edges 
by current XSPEC models. While these absorption mod- 
els are generally adequate for EPIC PN spectra, the data 
presented here has an exceptionally high signal-to-noise 
ratio so that very small discrepancies between data and 
model become apparent. We confirmed this interpreta- 
tion with the RGS spectra which clearly show a complex 
structure with several narrow absorption features near 
the K-shell absorption edges of oxygen and neon. This 
edge structure is not adequately described by any of the 
absorption models in XSPEC which model the edges as 
an exponential with a low-energy cut-off. The structure 
of the oxygen and neon K-edge s in the R GS spectra is 
similar to the structure found bv lJuett et a l. (2004, 200(|) 
for several X-ray binaries (including 4U 1636-536) and 
interpreted as absorption by neutral and ionized oxygen 
and neon in the interstellar medium. We find no emission 
lines in the RGS spectra that could explain the residu- 
als in the EPIC PN spectra. Because we are mainly 
interested in the spectrum at higher energies near the 
Fe K-a line, we simply excluded the data near the ab- 
sorption edges with the largest residuals. The following 
energy ranges were excluded from the fits: < 0.5 keV 
(O K-edge at 0.54 keV), 0.78-0.99 keV (Ne K-edge at 
0.87 keV), 1.13-1.45 keV (Mg K-edge at 1.3 keV), and 
1.75-1.87 keV (Si K-edge at 1.8 keV). 

In order to determine an appropriate model for the 
continuum emission, we initially fitted the spectra while 
excluding the 3-10 keV energy range. This excludes any 
contribution from a relativistically broadened, fluores- 
cent Fe K-a line which, accord ing to the diskline model 
in XSPEC (|Fabian et al.ll989f) . is limited to the 3-10 keV 
energy range for accretion disks around neutron stars. 
The X-ray spectra of LMXBs are typically characterized 
by a hard component interpreted as Comptonization of 
soft photons by a hot corona, a soft component thought 
to be blackbody radiation from the accretion disk, and a 
broad Fe K-a line at 6.4 keV (e.g. iBarret et alll2000f) . 
To fit the X-ray spectrum of the continuu m emission 
we th erefore combined a compTT component (!Titarchui3 
|1994|) for a Comptonizing corona or boundary layer, a 
disk bb disk blackbody component fe.g. lMakishima et al.l 
Il986f) . and a vphabs component for photoelectric absorp- 
tion with variable abundances (model 1). In the compTT 
model we used the approximation for a disk geometry 
and fixed the redshift at 0. In the vphabs model we only 
varied the elemental abundances for O, Ne, Si, and Fe. 
Abundances for the other elements were not well con- 
strained by the fit and w e re fix ed at their solar values 
according to iWilms et al.l ()2000D . The spectra from the 
three observations are reasonably well fitted by model 1 
with xl values of 1.44, 1.80, and 1.28 (Table[J). 
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TABLE 1 

Fits of Various XSPEC Models to the X-ray Spectra of 4U 1636-536 



Obs. 1 Obs. 2 Obs. 3 
Spectral Model xl (xVdof) xl (x^dof) xg (xVdof) 

1. vphabs*(compTT+diskbb) 1.44 (172.5/120) 1.80 (178.2/99) 1.28 (122.8/96) 

2. vphabs*(compTT+diskbb+bbodyrad) 1.36 (161.0/118) 1.43 (139.0/97) 0.97 (91.4/94) 

3. vphabs*(compTT+diskbb+bbodyrad+diskline) 1.23 (314.2/255) 1.22 (285.3/234) 1.18 (272.4/231) 

4. vphabs*(compTT+diskbb+bbodyrad+diskline+diskliiie) 1.23 (312.4/254) 1.21 (281.1/233) 1.19 (273.8/230) 

Note. — values for fits of various XSPEC models to the EPIC PN and HEXTE spectra of 4U 1636-536. Shown are the 
reduced (xi — X^/dof), the x^i and the number of degrees of freedom (dof). For models 1 and 2 the energy range 3-10 keV 
was excluded from the fit. 



We find that the fit can be significantly improved 
by adding a second blackbody component (bbodyrad) 
with a temperature of ~200 eV (model 2). The ad- 
ditional component improves the for the three ob- 
servations to 1.36, 1.43, and 0.97, respectively. A sec- 
ond blackbody compo nent with a similar temperature 
was also introduced bv lFiocchi et al.l (l2006l ) in their fit- 
ting of BeppoSAX spectra of 4U 1636-536 obtained dur- 
ing the high/soft state. The region associated with 
the blackbody component has an apparent radius of 
~100 km, which suggest that this emission originates 
from the accretion disk. We note that the diskbb model 
(|Makishima et all 11986^ was constructed for accretion 
disks around black holes and may not correctly describe 
the temperature profile in accretion disks around neu- 
tron stars which are subjected to significant central illu- 
mination from the stellar surface. The second blackbody 
component may represent deviations of the temperature 
profile in the neutron star accretion disk from the diskbb 
model. No further improvement in the fit was achieved 
by adding a third blackbody component or other com- 
mon spectral components. The remaining discrepancy 
between data and model is probably caused by small cal- 
ibration uncertainties that are noticeable because of the 
very high signal-to-noise ratio of the data and, as men- 
tioned above, by the overly simplified modeling of some 
absorption edges. 

For our further analysis we included the data in the 
3-10 keV energy range which were previously excluded 
when fitting models 1 and 2. The data in this energy 
range show a clear excess over the continuum flux pre- 
dicted by model 2 which is likely caused by relativistically 
broadened Fe K-a line emission. In order to fit the pro- 
file of this F e K-a line, we add ed to model 2 a diskline 
component (jFabian et al.lll989[ l which describes relativis- 
tically broadened line emission from a neutron star ac- 
cretion disk (model 3). All parameters of the diskline 
model were allowed to vary during the fit. We initially 
attempted to vary only the diskline parameters and fix 
all other parameters at their values from model 2, but we 
found that a better fit can be achieved by allowing the 
continuum model parameters to vary as well. This is an 
indication that the spectrum outside the 3-10 keV energy 
range is insufficient to correctly determine all continuum 
model parameters. Model 3 provides a good fit to the 
data with xt values of 1.23, 1.22, and 1.18, respectively. 

Figure[T]shows the observed spectra for the three obser- 
vations, the best fit with model 3, the individual model 
components, and the residuals between data and model. 
The excess of the line emission over the continuum model 
is shown in Figure [2l The figure clearly shows a broad 



and asymmetric Fe K-a line with a peak flux of 5-8% 
over the continuum. The line profile is well described 
by the diskline model. The best-fit parameters of the 
continuum and line components for model 3 are shown 
in Tables [5] and [3] Also shown are the uncertainties of 
the parameters at a 90% confidence level (95% for up- 
per and lower limits). When calculating the confidence 
regions, all continuum and line parameters were allowed 
to vary. We found that the commonly used error and 
steppar commands in XSPEC frequently failed to con- 
verge or underestimated the uncertainties. We therefore 
calculated the confidence regions for many of the param- 
eters by manually searching for the parameter value that 
produced the appropriate change in x^. 

The continuum parameters for model 3 are gener- 
ally consistent between the three observations, although 
the fit indicates a higher disk temperature and stronger 
Comptonization component for the later observations. 
We find a significant difference in the rest-frame energy 
Eq of the Fe K-a line. For observations 1 and 3, Eq is 
consistent with 6.4 keV for weakly ionized iron, whereas, 
for observation 2, Eq is close to 7.0 keV for highly ion- 
ized Fe XXVI. The power-law index /? of the emissivity 
profile is in the typically range found for iron lines from 
black hole accretion disks. The diskline model provides 
lower limits on the disk inclination of 81° for observations 
1 and 3 and 64° for observation 2. A lower limit of 81° 
is clearly inconsistent with the fact that 4U 1636-536 is 
not an eclipsing sys tem as well as with the results by 
ICasares et al.l ()2006l) who constrained the orbital incli- 
nation to 36-74° using phase-resolved spectroscopy. The 
high limit on the disk inclination for observations 1 and 
3 is likely the result of a slightly broader line profile com- 
pared to observation 2. Our fit with the diskline model 
suggests that the line profile is broader than physically 
possible for a single Fe K-a emission line. One possi- 
ble explanation for a broader line profile is a smaller 
radius of the inner disk edge Rm- However, the best- 
fit value of Rin for observations 1 and 3 is already at 
6Rg {Rg = GM/c^), the lower limit of the diskline 
model and the radius of the innermost stable circular or- 
bit (ISCO) for a non-rotating neutron star. An inner disk 
radius much smaller than QRg, while possible for rotat- 
ing black holes, is unlikely for the accretion disk around 
the neutron star in 4U 1636-536 which is rotating at a 
rate of 581 Hz. The most likely explanation for a broader 
line profile is the presence of several blended Fe K-a lines 
with different rest-frame energies. The presence of iron 
in more than one ionization state is already indicated 
by the difference in Eq between the three observations 
(6.30 keV and 6.43 keV vs. 7.06 keV). 
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Fig. 1.— EPIC PN and HEXTE spectra of 4U 1636-536 for the 
three observations. The upper panels show the observed flux {er- 
ror bars), the best fit with model 3 {solid line), and the individual 
additive model components {dotted lines). The individual model 
components are from loft to right: bbodyrad, diskbb, compTT, and 
diskline. The lower panels show the residuals between the ob- 
served flux and the model divided by the Icr error of each data 
point. 
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Fig. 2. — Relativistic iron line profiles in 4U 1636-536. Shown 
are the observed flux {error bars) and the model flux {solid line) 
normalized to the continuum flux from model 3 (i.e. the model flux 
without the diskline component). 

To test whether blended Fe K-a lines are a viable ex- 
planation, we added a second diskline component to 
the model (model 4). Because the iron lines are broad 
and significantly blended, the parameters of the two 
diskline components are strongly correlated, and it is 
not possible to fit all line parameters independently. We 
therefore fixed the line energy of the diskline compo- 
nents at 6.4 keV and 7.0 keV, respectively, and linked 
the disk inclination parameters. We also found that the 
outer disk radius Rout is poorly constrained by the fit, 
so we kept it fixed at 1000 Rg. The continuum parame- 
ters, which were allowed to vary simultaneously with the 
line parameters, did not change significantly compared 
to model 3. Adding the second line component does not 
significantly improve the fit (Table [T]), but it does result 
in more reasonable line parameters (Table [3]). The lower 
limits on the disk inclination change to 64° and 65°, re- 
spectively, which, unlike for model 3, are no w consistent 
with the 36-74° constraint by Casar es et all (|I006). The 
inner radius Rin for the 7.0 keV component appears to 
be generally smaller than for the 6.4 keV component as 
is expected for an accretion disk in which the tempera- 
ture increases towards the center. The fiux ratio of the 
7.0 keV and the 6.4 keV component increases from 0.9 
to 1.6 to 4.9 between the three observations, indicating 
that iron is generally in a higher ionization state for the 
later observations. This is consistent with the increase in 
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X-ray flux between observations. We conclude that the 
presence of iron at different ionization states is a viable 
explanation for the unusually broad line profiles and that 
at least two blended Fe K-a lines are necessary to obtain 
consistent line parameters. 

4. DISCUSSION 

We have analyzed X-ray spectra of the neutron star 
LMXB 4U 1636-536 obtained with XMM-Newton and 
RXTE in 2005, 2007, and 2008. The very high signal- 
to-noise ratio of the spectra allowed us to clearly de- 
tect a broad, relativistic Fe K-a line from the inner ac- 
cretion disk. The line is significantly broader than the 
asymmetric Fe K- g lines recently found in other neu- 
tron star LMXBs jBhattacharvva fc Strohmaveil l2007t 
ICackett et al.l 120081 ) . The broader fine profile is likely 
the result of a high disk inclination in 4U 1636-536. The 
inclination angles derived from iron line profiles in other 
neutron star LMX Bs have so far been co mparatively low. 
As pointed out by ICackett et al.l ()2008) , this is likely a 
selection effect because the narrower lines in low inclina- 
tion systems are more easily detectable. With the high 
signal-to-noise ratio of the 4U 1636-536 spectra, we have 
now been able to measure the broader line profile in a 
high inclination LMXB. Our analysis of the Fe K-a line 
profile places a lower limit of 64° on the disk inclination 
in 4U 1636-536. This limit is consistent with the 36- 
74° co nstraint on the orbital inclination bv lCasares et al.l 
()2006l ) and with the non-detection of eclipses. 

When fitting the iron line profile with a single rela- 
tivistic line component, we find a significant difference 
in the rest-frame line energy between the three obser- 
vations. This difference is likely caused by a change in 
the ionization profile of the disk related to the change 
in X-ray luminosity between the three observations. The 
line energy derived for the second observation is close to 
6.97 keV, the K-a line energy of Fe XXVI, which suggests 
that Fe XXVI contributes significantly to the observed 
line profile. According to our fit with a disk-blackbody 
model, the highest temperature in the disk is ~0.9 keV. 
Because plasma in thermal equilibrium at this tempera- 
ture does not contain a significant fraction of Fe XXVI, is 
it evident that the ionization profile in the disk is strongly 
affected by photoionization. 

We find that the Fe K-a line profile for two of the 
observations is too broad to be adequately described by 
a single relativistic emission line with physically reason- 
able values of disk inclination and inner disk radius. The 
broader than expected line profile can be explained by 
overlapping K-a lines from iron in different ionization 
states. The presence of multiple ionization states is also 
indicated by the difference in the fitted line energy be- 
tween the three observations. It is evident that multi- 
ple line components need to be considered to adequately 
model the relativistic iron line profiles in neutron star 
LMXBs. We obtained reasonable line parameters when 
fitting two iron lines with different rest-frame energies. 
This is obviously an oversimplification, since many ion- 
ization states are probably contributing to the iron line 
profile. However, because the lines are broad and overlap 
considerably, it is not possible to constrain the param- 
eters of all line components independently. Even with 
only two line components, the fit parameters are already 
strongly correlated, leading to larger parameter uncer- 



tainties than a fit with a single line component. In order 
to adequately model the contribution from multiple ion- 
ization states when fitting relativistic Fe K-a line pro- 
files, improved models are needed that can predict the 
ionization profile in the accretion disk. 

The Fe K-a line profiles in neutron star LMXBs can 
in principle be used to place upper limits on the neu- 
tron star radius by constraining the inner disk radius 
(jCackett et al.l[2008l ) . Previously, these line profiles have 
been fitted with only a single line component. However, 
if more than one ionization state contributes significantly 
to the Fe K-a emission, the line profile will be broader 
than for a single relativistic line, and a fit with a sin- 
gle line model may underestimate the inner disk radius 
and thus the limit on the neutron star radius. This can 
be clearly seen for two of our observations for which the 
upper limit on Rin increases from 6.3i?g to 11.9i?g and 
10.9i?g, respectively, when the iron line profile is fitted 
with two line components instead of a single line compo- 
nent. In contrast, the upper limit on Rin for the second 
observation decreases from 13.3i?g to 9.8Rg. It is evident 
that the constraints on the inner disk radius strongly de- 
pend on the assumptions made about the contributing 
ionization states of iron and that a better understanding 
of the ionization profile in the disk is needed to obtain re- 
liable limits on the neutron star radius. The contribution 
of multiple ionization states may also be important for 
the interpretation of some iron line profiles in accreting 
black holes. We note that the disk inclination and in- 
ner disk radius are strongly anti-correlated when fitting 
broad iron line profiles, which can lead to large uncer- 
tainties of the two parameters. Prior knowledge of the 
disk inclination can significantly reduce the uncertainty 
of the inner disk radius^ 

It was shown by iLaurent fc TitarchukI ()2007f ) that 
broad iron line features can also be produced by Comp- 
ton scattering of line photons in a strong outflow. The 
expected line proflles for this process are generally char- 
acterized by a narrow line at ~6.4-6.6 keV from fluores- 
cence in the outflow and a broad, redshifted component 
below 7 kcV from downscattering of the line photons. 
These line profiles differ qualitatively from those found in 
4U 1636-536 which show significant emission above 7 keV 
and no narrow line component (Figure [2|) . We also note 
that Compton scattering in an outflow only contributes 
significantly to the iron line emission if the optical depth 
is at least of order unity, which requires a mass outfiow 
rate on the order of the Eddington mass accretion rate. 
The X-ray luminosity we observed in 4U 1636-536 was 
only ~5% of the Eddington luminosity, which suggests 
that the rate of any outfiow in 4U 1636-536 was signifi- 
cantly below the Eddington mass accretion rate. It there- 
fore seems unlikely that a large fraction of the observed 
iron line emission was produced in an outflow. 
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TABLE 2 

Fit Parameters of the Continuum Components of Model 3 



T--'il T" a Wl O't" OT" 

IT Oil dlliU LUl 


Obs. 1 


Obs. 2 


Obs. 3 


vphabs:'' 








Nh (lO^i cm-2) 


3.63 ± 0.12 


3.83 ±0.10 


3.83 ±0.08 


O abundance 


1.37 ±0.04 


1.27 ± 0.03 


1.27 ± 0.03 


Ne abundance 


1.7± 0.2 


1.2 ±0.2 


1.3 ±0.2 


Si abundance 


1.7± 0.4 


2.0 ±0.4 


0.8 ±0.4 


Fe abundance 


1.49 ± 0.10 


1.47 ± 0.11 


1.48 ± 0.09 


compTT:'' 








To (keV) 


0.98 ± 0.02 


1.45 ± 0.02 


1.68 ± 0.03 


Tc (keV) 


14.2 ± 0.9 


> 9 


> 4 




2.34 ±0.12 


< 1.3 


< 2.7 


Normalization (10~^) 


9.4 ±0.7 


9±8 


25± 22 


Flux (lO"-*^" erg cm~^ 


6.1 


8.2 


10.0 


diskbb:'' 








T^sk (eV) 


640 ± 30 


780 ± 20 


900 ± 20 


Normalization 


150 ± 30 


200 ± 13 


152 ±7 


RdiskV cos i (km) 


7.3 ±0.7 


8.5 ±0.3 


7.4 ± 0.2 


Flux (10-1° erg cm'^ s^i) 


4.5 


13.5 


18.3 


bbodyrad:'' 








Tib (eV) 


174 ±9 


177 ± 7 


193 ± 6 


Normalization (10^) 


29 ± 6 


29 ± 8 


23 ±5 


Rbb (km) 


102 ± 11 


102 ± 14 


91 ± 10 


Flux (lO--"^" erg cm"^ s~^) 


1.9 


2.0 


2.3 


Note. — Uncertainties are given at 


a 90% confidence level and limits at a 95% confidence level. Fluxes 



for individual model components are unabsorbed and for the energy range 0.5 -10 keV. RdiikVcosi and R^t 
were calculated from the model normalizations assuming a distance of 6 kpc II Gallowav et al.ll2006l) . 
^ Photoelectric absorption model with variable abundances. Nh is the neutral hydrogen column density. 
Elemental abundances are given relative to their solar values according to Wilms et al. (2000). Abundances 
for elements not shown were fixed at their solar values. 

^ Comptonization model according to Titarchuk (1994). To is the temperature of the soft seed photons, 
Tq the temperature of the Comptonizing plasma, and tq the optical depth. A disk geometry was assumed 
for the model, and the redshift was fi xed at 0. 

Disk-blackbody model (e.g. jMakishima et anil986l). T^ia k is the temperature at the inner disk edge, 
Rdisk the apparent inner disk radius fsee lKubota et al.lll998l) . and i the disk inclination. 

Blackbody model. T^b is the blackbody temperature and R^b the apparent radius of the emitting region. 



TABLE 3 

Fit Parameters of the Iron Line Components of Models 3 and 4 



Obs. 1 Obs. 2 Obs. 3 

Model 3: 



Eo (keV) 


6.30 ±0.09 


7.06 ±0.10 


6.43 ± 0.09 


/9 


-2.72 ± 0.11 


-2.32 ± 0.22 


-2.73 ±0.13 


R,„ (GM/c^) 


< 6.3 




< 13.3 




< 6.3 


Rout {GM/c^) 


oon+2900 




2101180 


i n 


> 81 






> 64 




> 81 


Normalization (10 ^) 


1.49 ±0.13 


0.91 ± 0.20 


2.00 ± 0.20 


Equivalent width (eV) 


215 






98 




140 


Flux (IQ-l^ erg cm"^ s'^) 


15.7 






10.2 




21.5 


Model 4: 














Eo (keV) 


(6.4) 


(7.0) 


(6.4) 


(7.0) 


(6.4) 


(7.0) 


/9 


-2.7 ± 0.2 


< -2.7 


< -2.6 


9 Y+0.3 
^•'-1.6 


< -2.7 


-2.7 ± 0.2 


R,„ {GM/c^) 


< 17.8 


< 11.9 


q 0+5.6 


< 9.8 


< 33.6 


< 10.9 


Rout (GMjc^) 


(1000) 


(1000) 


(1000) 


(1000) 


(1000) 


(1000) 


i (°) 


> 64 


> 64 


> 65 


> 65 


> 64 


> 64 


Normalization (10 ^) 


Q 7+0.3 
'^•'-0.2 


6+°-2 
'-'■"-0.4 


5+°-'' 
'-'■^-0.3 




0.3 ±0.2 


1.3± 0.2 


Equivalent width (eV) 


108 


131 


54 


84 


23 


100 


Flux (IQ--^^ erg cm"^ s~^) 


7.7 


7.0 


5.7 


8.9 


3.0 


14.6 


Note. — Eq is the rest-frame 


energy of the 


iron emission line, /3 the power-law index of the emissivity 



dependence on radius, and Rout the inner and outer disk radius (in units of GM/c^ with M being the 
neutron star mass), and i the disk inclination. Also shown is the equivalent line width and the integrated 
flux for each diskline component. Uncertainties are given at a 90% confidence level and limits at a 95% 
confidence level. Values in parenthesis indicate that the parameter was fixed during the fit. 
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